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I. INTRODUCTION 

In recent years, considerable effort has been directed 
towards identifying and classifying phase transitions far 
from thermal equilibrium. Such nonequilibrium transi- 
tions can be found in a wide variety of problems in bi- 
ology, chemistry, and physics. Examples include pop- 
ulation dynamics, the spreading of epidemics, surface 
chemical reactions, catalysis, granular flow, traffic jams 
as well as growing surfaces and interfaces (see, e.g., 
[1, 2, 3, 4, 5, 6, 7, 8]). Nonequilibrium phase transi- 
tions are characterized by large scale fluctuations and 
collective behavior in space and time very similar to the 
behavior at equilibrium critical points. 

A particularly interesting situation arises when an 
equilibrium or nonequilibrium many-particle system is 
defined on a randomly diluted lattice. Then, two distinct 
types of fluctuations are combined, viz. the dynamical 
fluctuations of the many-particle system and the static 
geometric fluctuations due to lattice percolation [9]. In 
equilibrium systems, their interplay gives rise to novel 
universality classes for the thermal [10, 11, 12] and quan- 
tum [13, 14, 15, 16] phase transitions across the lattice 
percolation threshold. 

In this paper, we investigate the interplay between dy- 
namical fluctuations and geometric criticality in nonequi- 
librium many-particle systems. We focus on a particu- 
larly well-studied type of transitions, the so-called ab- 
sorbing state transitions, that separate active, fluctuating 
steady states from inactive (absorbing) states in which 
fluctuations cease completely. The generic universality 
class for absorbing state transitions is directed percola- 
tion (DP) [17]. It is conjectured [18, 19] to be valid for 
all absorbing state transitions with scalar order parame- 
ter and no extra symmetries or conservation laws. In the 
presence of symmetries and/or conservation laws, other 
universality classes can be realized, such as the DPn class 



in systems with n symmetric absorbing states [20] . 

For definitencss, we consider the contact process [21], 
a prototypical system in the DP universality class. We 
show that the contact process on a randomly site or bond 
diluted lattice has two different nonequilibrium phase 
transitions: (i) a generic disordered DP transition at 
weak dilutions (below the lattice percolation threshold) 
driven by the dynamic fluctuations of the contact pro- 
cess and (ii) the transition across the lattice percolation 
threshold driven by the geometric criticality of the lat- 
tice. The former transition has been investigated for a 
number of years [22, 23, 24, 25]; it has recently reat- 
tracted considerable attention because it is governed by 
an exotic infinite-randomness fixed point [26, 27, 28, 29]. 
In contrast, the latter transition has received much less 
attention. 

Here, we develop a theory for the nonequilibrium tran- 
sition across the lattice percolation threshold by combin- 
ing classical percolation theory with the properties of the 
supercritical contact process on a finite-size cluster. We 
show that the critical point is characterized by ultraslow 
activated (exponential) dynamical scaling and accompa- 
nied by strong Griffiths singularities. The scaling sce- 
nario is qualitatively similar to the generic disordered 
DP transition, but with different critical exponent val- 
ues. To confirm the universality of this exotic scenario, 
we also investigate the generalized contact process with n 
(symmetric) absorbing states [20] . This is a particularly 
interesting problem because the generic transition of the 
disordered generalized contact process does not appear 
to be of infinite-randomness type [26, 27]. 

The paper is organized as follows. In Sec. II, we in- 
troduce our models, the simple and generalized contact 
processes on a randomly diluted lattice. We also discuss 
the phase diagrams. In Sec. Ill we briefly summarize the 
results of classical percolation theory to the extent nec- 
essary for our purposes. Section IV contains the main 
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part of the paper, the theory of the nonequilibrium tran- 
sition across the lattice percolation threshold. Section V 
is devoted to the question of the generality of the arising 
scaling scenario. We conclude in Sec. VI. A short ac- 
count of part of this work has already been published in 
Rcf. [30]. 

II. SIMPLE AND GENERALIZED CONTACT 
PROCESSES ON DILUTED LATTICES 

A. Contact process 

The clean contact process [21] is a prototypical sys- 
tem in the DP universality class. It is defined on a d- 
dimensional hypercubic lattice. (We consider d > 2 since 
we will be interested in diluting the lattice.) Each lat- 
tice site r can be active (infected, state A) or inactive 
(healthy, state I) . During the time evolution of the con- 
tact process which is a continuous-time Markov process, 
each active site becomes inactive at a rate p ( "healing" ) 
while each inactive site becomes active at a rate Am/ (2d) 
where m is the number of active nearest neighbor sites 
("infection"). The infection rate A and the healing rate 
fi are external parameters. Their ratio controls the be- 
havior of the contact process. 

For A <C /i, healing dominates over infection, and the 
absorbing state without any active sites is the only steady 
state of the system (inactive phase) . For sufficiently large 
infection rate A, there is a steady state with a nonzero 
density of active sites (active phase). These two phases 
are separated by a nonequilibrium phase transition in 
the DP universality class at a critical value (A//i)° of the 
ratio of the infection and healing rates. 

The basic observable in the contact process is the av- 
erage density of active sites at time t, 

p(t) = ^sJ2(Mt)) (i) 

r 

where n r (t) — 1 if the site r is active at time t and 
n r (t) = if it is inactive. L is the linear system size, 
and (...) denotes the average over all realizations of the 
Markov process. The longtime limit of this density (i.e., 
the steady state density) 

Pstat = lim p(t) (2) 

t — >oo 

is the order parameter of the nonequilibrium phase tran- 
sition. 



B. Generalized contact process 

Following Hinrichsen [20], we now generalize the con- 
tact process by introducing n different inactive states I& 
with k = 1 . . . n (n = 1 corresponds to the simple contact 
process). Here, k is sometimes called the "color" label. 



The time evolution is again a continuous-time Markov 
process. The first two rates are equivalent to those of 
the simple contact process: An active site can decay into 
each of the inactive states I& with rate fi/n, and a site 
in any of the inactive states becomes active at a rate 
Am/ (2c?) with m the number of active nearest-neighbor 
sites. To introduce competition between the different in- 
active states, we define a third rate: If two neighboring 
sites are in different inactive states, each can become ac- 
tive with a rate a. This last rule prevents the boundaries 
between domains of different inactive states from stick- 
ing together infinitely. Instead they can separate, leaving 
active sites behind. 

The properties of the clean generalized contact process 
have been studied in some detail in the literature [20, 31]. 
If the boundary activation rate a vanishes, the behavior 
becomes identical to the simple contact process for all n. 
(This becomes obvious by simply dropping the color label 
and treating all inactive sites as identical.) For a > 0, the 
system becomes "more active" than the simple contact 
process, and the universality class changes. In one space 
dimension, a phase transition exists for n = 1 (in the DP 
universality class) and for n = 2 (in the .^-symmetric 
directed percolation (DP2) class which coincides with the 
the parity-conserving (PC) class in one dimension [5]). 
For n > 3 the system is always in the active phase, and 
no phase transition exists at finite values of A, p and a . 

The generalized contact process in higher space dimen- 
sions presumably behaves in an analogous fashion: There 
is a DP transition for n = 1 while the properties for n > 1 
arc different. For sufficiently large n, the system is always 
active [32]. 



C. Lattice dilution 

We now introduce quenched site dilution by randomly 
removing each lattice site with probability p. (Bond di- 
lution could be introduced analogously.) As long as the 
vacancy concentration p remains below the lattice per- 
colation threshold p c , the lattice consists of an infinite 
connected cluster of sites accompanied by a spectrum of 
finite-size clusters. In contrast, at dilutions above p c , the 
lattice consists of disconnected finite-size clusters only. 

Figure 1 schematically shows the resulting phase dia- 
grams of the nonequilibrium process as a function of the 
infection rate A and dilution p, keeping the healing rate 
/i and the boundary activation rate <r, if any, constant. 
Depending on the properties of the clean undiluted sys- 
tem, there are two qualitatively different cases. 

(a) If the undiluted system has a phase transition at 
a nonzero critical infection rate A° , the active phase sur- 
vives for all vacancy concentrations below the percolation 
threshold, p < p c . It even survives at the percolation 
threshold p c on the critical percolation cluster because it 
is connected, infinitely extended, and its fractal dimen- 
sion Df is larger than unity. The critical infection rate A c 
increases with increasing dilution p to compensate for the 
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FIG. 1: (Color online:) Schematic phase diagrams for the 
simple and generalized contact processes on a diluted lattice 
in dimensions d > 2 as a function of dilution p and inverse 
infection rate A -1 (healing and boundary activation rates fj, 
and a are fixed). Case (a) applies to systems that display a 
phase transition at A° in the absence of dilution. There is a 
multicritical point (MCP) at (p c ,A,) separating the generic 
transition from the lattice percolation transition. Case (b) is 
for systems that are always active in the absence of dilution. 



summarize key results of percolation theory. 



III. CLASSICAL PERCOLATION THEORY 

Consider a regular lattice in d dimensions. If each lat- 
tice site is removed with probability p [34], an obvious 
question is whether or not the lattice is still connected 
in the sense that there is a cluster of connected (nearest 
neighbor) sites that spans the entire system. This ques- 
tion defines the percolation problem (see Ref. [9] for an 
introduction) . 

In the thermodynamic limit of infinite system volume, 
there is a sharp boundary between the cases of a con- 
nected or disconnected lattice. If the vacancy concen- 
tration p stays below the percolation threshold p Cl an 
infinite cluster of connected sites exists (with a proba- 
bility of unity). For p > p c , an infinite cluster does not 
exist, instead, the lattice consists of many disconnected 
finite-size clusters. 

The behavior of the lattice for vacancy concentrations 
close to the percolation threshold can be understood as 
a (geometric) continuous phase transition or critical phe- 
nomenon. The order parameter is the probability 
of a site to belong to the infinite connected percolation 
cluster. It is obviously zero in the disconnected phase 
(p > p c ) and nonzero in the percolating phase (p < p c ). 
Close to p c it varies as 



\P~P, 



I /3c 



(P < Pc) 



(3) 



missing neighbors, reaching A* at p c . The active phase 
cannot exist for p > p c because the lattice consists of 
finite-size clusters only, and the noncquilibrium process 
will eventually end up in one of the absorbing states on 
any finite-size cluster. Thus, in case (a), our system fea- 
tures two noncquilibrium phase transitions, (i) a generic 
(disordered) transition for dilutions p < p c , driven by the 
dynamic fluctuations of the noncquilibrium process and 
(ii) the transition across the lattice percolation threshold 
driven by the geometric criticality of the lattice. They 
arc separated by a multicritical point at (p c ,A„) which 
was studied numerically in Ref. [33] . 

(b) If the undiluted system is always active (as for the 
generalized contact process with a sufficiently high num- 
ber of inactive states), the phase diagram is simpler. The 
active phase covers the entire region p < p c for all A > 
(A* is formally zero) while the inactive phase exists in the 
region p > p c . There is no generic (disordered) nonequi- 
librium phase transition, only the transition across the 
lattice percolation threshold. 

The focus of the present paper is the nonequilibrium 
phase transition across the lattice percolation threshold 
that exists in both cases. In order to develop a theory for 
this transition, we combine classical percolation theory 
with the properties of the nonequilibrium process on a 
finite-size cluster. In the next section we therefore briefly 



where (3 C is the order parameter critical exponent of clas- 
sical percolation. Note that we use a subscript c to distin- 
guish quantities associated with the classical lattice per- 
colation problem from those of the noncquilibrium phase 
transitions discussed later. In addition to the infinite 
cluster, we also need to characterize the finite clusters 
on both sides of the transition. Their typical size, the 
correlation or connectedness length £ c diverges as 

i*~\p-Pc\- Va (4) 

with v c the correlation length exponent. The average 
mass S c (number of sites) of a finite cluster diverges with 
the susceptibility exponent j c according to 

Sc-\p-p c \-^. (5) 

The complete information about the percolation criti- 
cal behavior is contained in the cluster size distribution 
n s , i.e., the number of clusters with s sites excluding the 
infinite cluster (normalized by the total number of lat- 
tice sites). Close to the percolation threshold, it obeys 
the scaling form 

n s (A) = s- T *f(As°°). (6) 

Here A = p — p c , and r c and <r c are critical exponents. 
The scaling function f(x) is analytic for small x and has 



a single maximum at some x max > 0. For large it 
drops off rapidly 



cxp 
cxp 



(x > 0), 

l-l/d' 



(7) 

(x < 0), (8) 



where B\ and B2 are constants of order unity. All 
classical percolation exponents are determined by t c 
and <7 C including the correlation lengths exponent v c = 
(r c — l)/(da c ), the order parameter exponent (3 C = (r c — 
2)/cr c , and the susceptibility exponent j c = (3 — t c )/ct c . 

Right at the percolation threshold, the cluster size dis- 
tribution does not contain a characteristic scale. The 
structure of the critical percolation cluster is thus frac- 
tal with the fractal dimension being given by Df = 
d/(r c - 1). 



IV. NONEQUILIBRIUM TRANSITION ACROSS 
THE LATTICE PERCOLATION THRESHOLD 

A. Single-cluster dynamics 

To develop a theory of the nonequilibrium phase tran- 
sition across the lattice percolation threshold, we first 
study the nonequilibrium process on a single connected 
finite-size cluster of s sites. For definitcness, this section 
focuses on the simple contact process. The generalized 
contact process will be considered in Sec. V. 

The crucial observation is that on the percolation tran- 
sition line (for A > A*), the contact process is supercrit- 
ical, i.e., the cluster is locally in the active phase. The 
time evolution of such a cluster, starting from a fully ac- 
tive lattice, therefore proceeds in two stages: Initially, 
the density p s of active sites decays rapidly towards a 
metastable state (which corresponds to the steady state 
of the equivalent infinite system) with a nonzero density 
of active sites and islands of the inactive phase of linear 
size (see Fig. 2). This metastable state can then decay 
into the inactive (absorbing) state only via a rare collec- 
tive fluctuation involving all sites of the cluster. We thus 
expect the long-time decay of the density to be of ex- 
ponential form (suppressing subleading pre-exponcntial 
factors) , 



p s (t) ~ exp[-t/t s (s)) , 



(9) 



with a long lifetime t s that increases exponentially with 
the cluster size s 



t a (s) = i expL4(A)s] 



(10) 



for sufficiently large s. Here, to is some microscopic time 
scale. 

The lifetime increases the faster with s the further the 
cluster is in the active phase. This means, the prefactor 
A(X) which plays the role of an inverse correlation volume 




FIG. 2: (Color online:) Schematic of the metastable state 
of the supercritical contact process on a single percolation 
cluster. A and I denote active and inactive sites, and £J is 
the connected correlation length of the density fluctuations 
on the cluster. 



vanishes at the multicritical value A* and monotonically 
increases with increasing A. Close to the multicritical 
point, the behavior of A(X) can be inferred from scaling. 
Since A(A) has the dimension of an inverse volume, it 
varies as 



A{\) ~ (A - \*) v ' s 



111) 



where is the correlation length exponent of the multi- 
critical point and Df is the (fractal) space dimensionality 
of the underlying cluster. 

Note that (10) establishes an exponential relation be- 
tween length and time scales at the transition. Because 
the number of sites s of a percolation cluster is related 
to its linear size R s via s ~ R^ f , eq. (10) implies 



Int. ~i2 



D f 



(12) 



Thus, the dynamical scaling is activated rather than 
power-law with the tunneling exponent being identical 
to the fractal dimension of the critical percolation clus- 
ter, ip = Df. 

To confirm the above phenomenological arguments, we 
have performed extensive Monte-Carlo simulations of the 
contact process on finite-size clusters using clean one- 
dimensional and two-dimensional systems as well as di- 
luted lattices. Our simulation method is based on the al- 
gorithm by Dickman [35] and described in detail in Rcfs. 
[28, 29]. 

A characteristic set of results is shown in Fig. 3. It 
shows the time evolution of the contact process on sev- 
eral one-dimensional clusters of different size s, starting 
from a fully active lattice. The infection rate A = 3.8 
(we set fx = 1) puts the clusters (locally) in the ordered 
phase, i.e., it is supercritical, since the critical value in 
one dimension is A c = 3.298. All data are averages over 
10 5 independent trials. The double-logarithmic plot of 
density p s vs. time t in Fig. 3a clearly shows the two- 
stage time evolution consisting of a rapid initial decay 
(independent of cluster size) towards a metastable state 
followed by a long-time decay towards the absorbing state 
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FIG. 3: (Color online:) Contact process on one-dimensional 
clusters of size s, starting from a fully active lattice at A = 
3.8, \i = 1 which is in the active phase, (a) Double-logarithmic 
plot of density vs. time showing the two-stage time-evolution 
via a metastable state, (b) Log-linear plot demonstrating that 
the long-time decay is exponential. All data are averages over 
10 independent runs. 



which becomes slower with increasing cluster size. Re- 
plotting the data in log-linear form in Fig. 3b confirms 
that the long-time decay is exponential, as predicted in 

(9) . 

The lifetime t s of the contact process on the cluster can 
be determined by fitting the asymptotic part of the p s (t) 
curve to (9). Figure 4 shows the lifetime as a function 
of cluster size s for four different values of the infection 
rate A. Clearly, for sufficiently large clusters, the lifetime 
depends exponentially on the cluster size, as predicted in 

(10) . (The data for A = 3.4 which is very close to the 
bulk critical point of A c = 3.298 have not fully reached 
the asymptotic regime as can be seen from the remaining 
slight curvature of the plot.) By fitting the large-s be- 
havior of the lifetime curves to the exponential law (10), 
we obtain an estimate of the inverse correlation volume 
A. The inset of Fig. 4 shows this correlation volume as 
a function of the distance from the bulk critical point. 
In accordance with (11) it behaves as a power law. The 
exponent value of approximately 0.95 is in reasonable 



FIG. 4: (Color online:) Lifetime t s as a function of cluster 
size s for different values of the infection rate A. The other 
parameters are as in Fig. 3. The dashed lines are fits of the 
large-s behavior to the exponential dependence (10). Inset: 
Correlation volume A -1 as a function of the distance from 
bulk criticality. The dashed line is a power-law fit. 

agreement with the prediction v = 1.097 for our one- 
dimensional clusters. 

We have performed analogous simulations for various 
sets of two-dimensional clusters as well as finite-size di- 
luted lattices. In all cases, the Monte-Carlo results con- 
firm the phenomenological theory summarized in eqs. (9), 
(10), and (11). 

B. Steady-state density and density decay 

We now consider the full problem, the contact process 
on a diluted lattice close to the percolation threshold. 
To obtain observables of the entire system, we must sum 
over all percolation clusters. 

Let us start by analyzing static quantities such as the 
steady state density p s t of active sites (the order parame- 
ter of the nonequilibrium transition) and the spatial cor- 
relation length £j_ . Finite-size percolation clusters do not 
contribute to the steady-state density because the con- 
tact process eventually decays into the absorbing inactive 
state on any finite-size cluster. A steady-state density 
can thus exist only on the infinite percolation cluster for 
p < p c . For A > A*, the infinite cluster is supercritical, 
i.e., a finite fraction of its sites is active. Thus, the to- 
tal steady-state density is proportional to the number of 
sites in the infinite cluster, 

n ~P fVl ~ I \P-Pcf" (P<Pc) (iTl 

ftt-'-W-J o (p> Pc ) ■ 

Consequently, the order parameter exponent [3 of the 
nonequilibrium transition is identical to the correspond- 
ing exponent j3 c of the lattice percolation problem. 

The (average) spatial correlation length £j_ of the 
nonequilibrium process can be found using a similar ar- 
gument. On the one hand, the spatial correlations of the 
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contact process cannot extend beyond the connectedness 
length £ c of the underlying diluted lattice because differ- 
ent percolation clusters are completely decoupled. This 
implies £j_ < £ c . On the other hand, for A > A*, all sites 
on the same percolation cluster are strongly correlated in 
space, implying £_l > £ c - We therefore conclude 



clusters of size s = -A~ 1 ]n[B 1 A 1 / <,c t /(At)]. Insert- 
ing this into the integrand results in a power-law density 
decay 



p(t,A)~(t/t y d /*' (p> Pc ) 



(20) 



(14) 



and the correlation length exponent v± is also identical 
to its lattice percolation counterpart v c . 

We now turn to the dynamics of the nonequilibrium 
transition across the percolation threshold. In order to 
find the time evolution of the total density of active sites 
(starting from a completely active lattice), we sum over 
all percolation clusters by combining the cluster size dis- 
tribution (6) with the single-cluster time evolution (9). 
The total density is thus given by 

p{t,A) = J dssn s (A) Ps (t) 

dssn s (A) exp[— t/t s (s)] (15) 



The nonuniversal exponent z' is given by z' = 
(Ad/ ' B\)A~ 1 /' Tc ~ €± f i i- c -i it diverges at the critical 
point p = p c . 

In the percolating, active phase (p < p c ), the infinite 
percolation cluster contributes a nonzero steady state 
density p st (A) given by (13). However, the long-time 
approach of the density towards this value is determined 
by the slow decay of the metastable states of large finite- 
size percolation clusters. To estimate their contribution, 
we must use the expression (8) for the scaling function of 
the cluster size distribution. The resulting integral now 
reads 



p(t,A)-p st (A) 



ds t 



exp 



(t/t e As )] 



-(B 2 s\A\ 1/ ^) 1 - 1/d 
(21) 



In the following, we evaluate this integral at the transi- 
tion as well as in the active and inactive phases. 

Right at the percolation threshold, the scaling function 
in the cluster size distribution (6) is a constant, /(0), and 
(15) simplifies to 



We again apply the saddle-point method to find the lead- 
ing low-time behavior of this integral. Minimizing the 
exponent shows the main contribution coming from clus- 
ters of size s = -A' 1 ln^lAI 1 /^^ - l)/(Atd)}. By 
inserting this into the integrand, we find a noncxponen- 
tial density decay of the form 



p(t,0)~ I ds s 1 -^ exp[-(t/t e As )} . (16) p(t, A) - p st (A) ~ e -[W z ") H*/*o)] 



To estimate this integral, we note that only suffi- 
ciently large clusters, with a minimum size of s m i n (t) = 
A^ 1 ln(i/io), contribute to the total density at time t, 



ds s 1 " 



(17) 



The leading long-time dependence of the total density 
right at the percolation threshold thus takes the unusual 
logarithmic form 



p(t,0) ~ [ln(i/i )] 



-8 



(18) 



again reflecting the activated dynamical scaling, with the 
critical exponent given by 8 = t c — 2 = f3 c / (v c D f). 

In the disconnected, inactive phase (p > p c ) we need 
to use expression (7) for the scaling function of the clus- 
ter size distribution. The resulting integral for the time 
evolution of the density reads 



p{t,A) ~ J dss 1 ^" exp[-B lS A 1/a " - (t/t e As )}. 



(19) 



For long times, the leading behavior of the integral can 
be calculated using the saddle-point method. Minimiz- 
ing the exponent of the integrand shows that the main 
contribution at time t to the integral (19) comes from 



(P < Pc) ■ 

(22) 

Here, z" = (Ad/B 2 )\A\~ 1 ^ a ~ ^ is another nonuni- 
versal exponent which diverges at the critical point. 

The slow noncxponcntial relaxation of the total den- 
sity on both sides of the actual transition as given in (20) 
and (22) is characteristic of a Griffiths phase [36] in the 
contact process [37]. It is brought about by the compe- 
tition between the exponentially decreasing probability 
for finding a large percolation cluster off criticality and 
the exponentially increasing lifetime of such a cluster. 
Note that time t and spatial correlation length £j_ enter 
the off-critical decay laws (20) and (22) in terms of the 
combination ln(i/io)/£j/ again reflecting the activated 
character of the dynamical scaling. 



C. Spreading from a single seed 

After having discussed the time evolution of the den- 
sity starting from a completely infected lattice, we now 
consider the survival probability P s (t) for runs starting 
from a single random seed site. To estimate P s (t), we 
note that the probability of a random seed site to belong 
to a cluster of size s is given by sn s (A). The activity of 
the contact process is confined to this seed cluster. Fol- 
lowing the arguments leading to (9) , the probability that 
this cluster survives is proportional to exp(— t/t s ). The 
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average survival probability at time t can thus be written 
as a sum over all possible seed clusters, 



P„(t,A) ~ J dssn s (A) cxp[- 



■t/t s (s)} 



(23) 



This is exactly the same integral as the one governing the 
density decay (15). We conclude that the time depen- 
dence of the survival probability for runs starting from a 
single seed is identical to the time evolution of the density 
when starting from a fully infected lattice, as is expected 
for the contact process under very general conditions (see, 
e.g., Rcf. [5]). 

To determine the (average) total number N(t) of ac- 
tive sites in a cloud spreading from a single seed, we 
observe that a supercritical cloud initially grows ballis- 
tically. This means its radius grows linearly with time, 
and the number of active sites follows a power law. This 
ballistic growth stops when the number of active sites is 
of the order of the cluster size s. After that, the number 
of active sites stays approximately constant. The num- 
ber N a (t) of active sites on a percolation cluster of size s 
is thus given by 



N s (t) 



(t/to) 1 



(t < ti(a)) 
(t > ti(s)) 



(24) 



where U(s) ~ R s (s) ~ toS 1 / D f is the saturation time of 
this cluster. Note that N s decays to zero only after the 
much longer cluster lifetime t s {s) = ioexpL4(A)s] given 
in (10). 

We now average over all possible positions of the seed 
site as in (23). This yields 



N{t, A) ' 

with s m j n ~ A^ 1 ln(t/t c 
easily evaluated, giving 

N(t,0) ~ t D f® 



dssn s (A) N s (t) 



(25) 



At criticality, this integral is 

) = ft./* _ (2 6 ) 



The lower bound of the integral (i.e., the logarithmically 
slow long-time decay of the clusters) produces a sub- 
leading correction only. Consequently, we arrive at the 
somewhat surprising conclusion that the initial spread- 
ing follows a power-law and is thus much faster than 
the long-time density decay. In contrast, at the infinite- 
randomness critical point governing the generic (p < p c ) 
transition, both the initial spreading and the long-time 
decay follow logarithmic laws [26, 27, 28, 29]. Note 
that a similar situation occurs at the percolation quan- 
tum phase transition in the diluted transverse-field Ising 
model [13] where the temperature-dependence of the cor- 
relation length does not follow the naively expected log- 
arithmic law. 



D. External source field 

In this subsection we discuss the effects of spontaneous 
activity creation on our noncquilibrium phase transition. 



Specifically, in addition to healing and infection, we now 
consider a third process by which an inactive site can 
spontaneously turn into an active site at rate h. This 
rate plays the role of an external "source field" conjugate 
to the order parameter. 

To find the steady state density in the presence of 
such a source field, we first consider a single percola- 
tion cluster. As before, we are interested in the super- 
critical regime A > A*. At any given time t, a clus- 
ter of size s will be active (on average), if at least one 
of the s sites has spontaneously become active within 
one lifetime t s (s) = toe As before t, i.e., in the interval 
[t — t s (s),t]. For a small external field h, the average 
number of active sites created on a cluster of size s is 
M s (h) = hst s (s) = hst$e As . This linear response expres- 
sion is valid as long as M s -C s. The probability w s (h) 
for a cluster of size s to be active in the steady state is 
thus given by 



w s (h) 



M 8 (h) 
1 



(M s (h) < 1) 
(M s (h) > 1) 



(27) 



Turning to the full lattice, the total steady state den- 
sity is obtained by summing over all clusters 



Prt (ft, A) ~ / dssn s (A)mm[l,M s (h)] 



(28) 



This integral can be evaluated along the same lines as 
the corresponding integral (15) for the time-evolution of 
the zero-field density. For small fields h, we obtain 



MM) 

Art (ft, A) 
<Wft,A) 



[Hh /h)}-~ 5 

(ft/fto) d/ *' 

[(d/z")ln(h/h )] 1 - 1/d 



(p = p c ), (29) 
(p > p c ), (30) 

(p < p c ) , (31) 



where Sp st (h, A) = p st (h, A)— p st (0, A) is the excess den- 
sity due to the field in the active phase and hg = 1/io- At 
criticality, p = p c , the relation between density Art an d 
field h is logarithmic because the field represents a rate 
(inverse time) and the dynamical scaling is activated. Off 
criticality, we find strong Griffiths singularities analogous 
to those in the time-dependence of the density. The ex- 
ponents z' and z" take the same values as calculated after 
eqs. (20) and (22), respectively. 



E. Scaling theory 

In Sections IV B and IV D, we have determined the 
critical behavior of the density of active sites by explicitly 
averaging the single cluster dynamics over all percolation 
clusters. The same results can also be obtained from 
writing down a general scaling theory of the density for 
the case of activated dynamical scaling [28, 29]. 

According to (13), in the active phase, the density is 
proportional to the number of sites in the infinite per- 
colation cluster. Its scale dimension must therefore be 



identical to the scale dimension of which is j3 c /v c . 
Time must enter the theory via the scaling combination 
ln(t/to)b^ with the tunneling exponent given by ip = Df 
and b an arbitrary length scale factor. This scaling com- 
bination reflects the activated dynamical scaling, i.e., the 
exponential relation (12) between length and time scales. 
Finally, the source field h, being a rate, scales like inverse 
time. This leads to the following scaling theory of the 
density, 

p[A,\n(t/t )M(h /h)} = 

= bP"/"- plAb- 1 ^- Mit/toP'^Hho/h^} (32) 

This scaling theory is compatible with all our explicit 
results which can be rederived by setting the arbitrary 
scale factor b to the appropriate values. 



V. GENERALITY OF THE ACTIVATED 
SCALING SCENARIO 

In Section IV, we have developed a theory for the 
nonequilibrium phase transition of the simple contact 
process across the lattice percolation threshold and found 
it to be characterized by unconventional activated dy- 
namical scaling. In the present section, we investigate 
how general this exotic behavior is for absorbing state 
transitions by considering the generalized contact pro- 
cess with several absorbing states. 

This is a particularly interesting question because the 
generic transitions (p < p c ) of the diluted simple and gen- 
eralized contact processes appear to behave differently 
The generic transition in the simple contact process has 
been shown to be of infinite-randomness type with ac- 
tivated dynamical scaling using both a strong-disorder 
renormalization group [26, 27] and Monte-Carlo simu- 
lations [28, 29]. In contrast, the strong-disorder renor- 
malization group treatment of the disordered generalized 
contact process [27] suggests more conventional behavior, 
even though the ultimate fate of the transition could not 
be determined. 

To address the same question for our transition across 
the lattice percolation threshold, we note that any dif- 
ference between the simple and the generalized contact 
processes must stem from the single-cluster dynamics be- 
cause the underlying lattice is identical. In the following 
we therefore first give heuristic arguments for the single- 
cluster dynamics of the supercritical generalized contact 
process and then verify them by Monte-Carlo simula- 
tions. 

If the percolation cluster is locally in the active phase 
(A > A*), the density time evolution, starting from a fully 
active lattice, proceeds in two stages, analogously to the 
simple contact process. There is a rapid initial decay to 
a metastable state with a nonzero density of active sites 
and finite-size islands of each of the inactive phases (see 
Fig. 5) . For this metastable state to decay into one of the 
n absorbing configurations, all sites must go into the same 
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FIG. 5: (Color online:) Schematic of the metastable state of 
the supercritical generalized contact process with two inactive 
states on a single percolation cluster. A denotes the active 
state, and Ii and I2 are the inactive states. £J is the connected 
correlation length of the density fluctuations on the cluster. 

inactive state which requires a rare large density fluctua- 
tion. Let us assume for definiteness that the decay is into 
the Ii state. The main difference to the simple contact 
process considered in Sec. IV A is that sites that are in in- 
active states I2 . . . I„ cannot directly decay into Ii . This 
means, each of the inactive islands in states I2 . . . I n first 
needs to be "eaten" by the active regions before the entire 
cluster can decay into the Ii state. This can only happen 
via infection from the boundary of the inactive island and 
is thus a slow process. However, since the characteristic 
size of the inactive islands in the metastable state is finite 
(it is given by the connected density correlation length 

on the cluster), this process happens with a nonzero 
rate that is independent of the size s of the underlying 
percolation cluster (for sufficiently large s). 

The decay of the metastable state into one of the ab- 
sorbing states is therefore brought about by the rare col- 
lective decay of a large number of independent correla- 
tion volumes just as in the simple contact process. As 
a result, the lifetime t s (s) depends exponentially on the 
number of involved correlation volumes, i.e., it depends 
exponentially on the cluster size s. We thus find that 
the long-time density decay of the generalized contact 
process on a single large percolation cluster is governed 
by the same equations (9) and (10) as the decay of the 
simple contact process. 

To verify these phenomenological arguments, we have 
performed large-scale Monte-Carlo simulations of the 
generalized contact process with two and three absorb- 
ing states on clean and disordered one-dimensional and 
two-dimensional lattices. In all cases, we have first per- 
formed bulk simulations (spreading from a single seed) to 
find the bulk critical point. An example is shown in Fig. 
6, details of the bulk critical behavior will be reported 
elsewhere. 

After having determined the critical point, if any, we 
have selected several parameter sets in the bulk active 
phase and studied the long-time density decay of the 
generalized contact process on finite size clusters. As ex- 
pected, the decay proceeds via the two stages discussed 
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FIG. 6: (Color online:) Bulk phase transition of the gener- 
alized contact process with two absorbing states in d = 1 
measured via spreading from a single seed: Number N of ac- 
tive sites vs. time t for different healing rates /i. The infection 
and boundary activation rates are fixed, A = a — 1, and the 
data are averages over 10 6 runs. The critical point appears 
to be close to fi — 0.628 in agreement with [20]. 



above. As in Sec. IV A, we extract the lifetime t s from 
the slow exponential long-time part of the decay. Two 
characteristic sets of results are shown in Fig. 7. The fig- 
ure confirms that the lifetime of the generalized contact 
process on a finite-size cluster depends exponentially on 
the number of sites in the cluster, as given in (10). We 
have obtained analogous results for all cases investigated, 
verifying the phenomcnological theory given above. 

Because the long-time dynamics of the generalized con- 
tact process on a single supercritical cluster follows the 
same behavior (9) and (10) as that of the simple contact 
process, we conclude that its noncquilibrium transition 
across the percolation threshold will also be governed by 
the theory developed on Sec. IV. In other words, the lat- 
tice percolation transitions of the simple and generalized 
contact processes belong to the same universality class, 
irrespective of the number n of absorbing states. 



VI. CONCLUSIONS 

In this final section of the paper, we first summarize 
our results, discuss their generality, and relate them to 
the behavior of certain quantum phase transitions on 
diluted lattices. We then compare the recently found 
infinite-randomness critical point at the generic transi- 
tion (j> < p c ) to the behavior at our lattice percolation 
transition. Finally, we relate our findings to a general 
classification of phase transitions with quenched spatial 
disorder [38]. 

To summarize, we have investigated absorbing state 
phase transitions on randomly diluted lattices, taking 
the simple and generalized contact processes as examples. 
We have focused on the noncquilibrium phase transition 
across the lattice percolation threshold and shown that it 
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FIG. 7: (Color online:) Lifetime t s as a function of cluster 
size s for the generalized contact process with two inactive 
states at different values of the healing rate fi. The infection 
and boundary activation rates are fixed, A = a = 1, and the 
data are averages over 10 6 runs, (a) d = 1 where the bulk 
system has a transition, see Fig. 6. (b) d = 2, where we do 
not find a bulk transition because the system is always active 
[32]. The dashed lines are fits of the large-s behaviors to the 
exponential law (10). 



can be understood by combining the time evolution of the 
supercritical nonequilibrium process on a finite-size clus- 
ter with results from classical lattice percolation theory. 
The interplay between geometric criticality and dynamic 
fluctuations at this transition leads to a novel univer- 
sality class. It is characterized by ultraslow activated 
(i.e., exponential) rather than power-law dynamical scal- 
ing and accompanied by a noncxponcntial decay in the 
Griffiths regions. All critical exponents of the nonequi- 
librium phase transition can be expressed in terms of the 
classical lattice percolation exponents. Their values are 
known exactly in two space dimensions and with good 
numerical accuracy in three space dimensions; they are 
summarized in Tabic I. Thus, our transition in d = 2 pro- 
vides one of the few examples of a noncquilibrium phase 
transition with exactly known critical exponents. 

The logarithmically slow dynamics (18), (29) at crit- 
icality together with the small value of the exponent 5 
make a numerical verification of our theory by Simula- 
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TABLE I: Critical exponents of the nonequilibrium phase 
transition across the percolation threshold in two and three 
space dimensions. 



Exponent 


d = 2 


d = 3 


P=0c 


5/36 


0.417 


v — V c 


4/3 


0.875 


i\> = Df = d — Pc/vc 


91/48 


2.523 


5 = P c /{v c D f ) 


5/91 


0.188 



tions of the full diluted lattice a very costly proposition. 
The results of recent Monte-Carlo simulations in two di- 
mensions [29] at p = p c are compatible with our theory 
but not yet sufficient to be considered a quantitative ver- 
ification. This remains a task for the future. 

The unconventional critical behavior of our nonequi- 
librium phase transition at p = p c is the direct result 
of combining the power-law spectrum (6) of cluster sizes 
with the exponential relation (12) between length and 
time scales. We therefore expect other equilibrium or 
nonequilibrium systems that share these two character- 
istics to display similar critical behavior at the lattice 
percolation transition. One prototypical example is the 
transvcrsc-ficld Ising model on a diluted lattice. In this 
system, the quantum-mechanical energy gap (which rep- 
resents an inverse time) of a cluster decreases exponen- 
tially with the cluster size. Consequently, the critical be- 
havior of the diluted transverse-field Ising model across 
the lattice percolation threshold is very similar to the 
one found in this paper [13]. Other candidates are mag- 
netic quantum phase transitions in metallic systems or 
certain superconductor-metal quantum phase transitions 
[39, 40, 41, 42], even though a pure percolation scenario 
may be hard to realize in metallic systems. 

Our work has focused on the nonequilibrium phase 
transition across the lattice percolation threshold. It 
is instructive to compare its critical behavior to that 
of the generic transition occurring for p < p c (see Fig. 
1). Hooyberghs et al. [26, 27] applied a strong disor- 
der renormalization group to the one-dimensional disor- 
dered contact process. They found an exotic infinite- 
randomness critical point in the universality class of the 
random-transverse field Ising model (which likely governs 
the transition for any disorder strength [43]). The same 
analogy is expected to hold in two space dimensions. Re- 
cently, these predictions were confirmed by large scale 
Monte-Carlo simulations [28, 29]. Our nonequilibrium 
transition across the lattice percolation threshold shares 
some characteristics with these infinite-randomness criti- 
cal points, in particular, the activated dynamical scaling 
which leads to a logarithmically slow density decay at 
criticality. 



However, the generic and percolation transitions are in 
different universality classes with different critical expo- 
nent values. Moreover, the initial spreading from a single 
seed is qualitatively different (logarithmically slow at the 
generic infinite-randomness critical point but of power- 
law type at our percolation transition). Finally, at the 
percolation transition the simple and generalized contact 
processes are in the same universality class while this 
does not seem to be the case for the generic transition 
[27]. 

The results of this paper are in agreement with a recent 
general classification of phase transitions with quenched 
spatial disorder and short-range interactions [38, 39]. It is 
based on the effective dimensionality d e s of the droplets 
or clusters. Three classes need to be distinguished: (a) If 
the clusters are below the lower critical dimension of the 
problem, d ff < d~ , the critical behavior is conventional 
(power-law scaling and exponentially weak Griffiths ef- 
fects). This is the case for most classical equilibrium 
transitions, (b) If d c ft = d~ , the dynamical scaling is ac- 
tivated and accompanied by strong Griffiths effects. This 
case is realized at the nonequilibrium transition consid- 
ered here as well as the generic transition of the disor- 
dered contact process. It also applies to various quan- 
tum phase transitions [13, 40, 44]. (c) If d e g > d~ , a 
single supercritical cluster can undergo the phase transi- 
tion independently of the bulk system. This leads to the 
smearing of the global phase transition; it occurs, e.g., in 
dissipative quantum magnets [45, 46] or in the contact 
process with extended defects [47]. 

In conclusion, our work demonstrates that absorbing 
state transitions on percolating lattices display unusual 
behavior. Interestingly, experimental verifications of the 
theoretically predicted critical behavior at (clean) ab- 
sorbing state transitions are extremely rare [48]. For in- 
stance, to the best of our knowledge, the only complete 
verification of directed percolation scaling was found very 
recently in the transition between two turbulent states 
in a liquid crystal [49]. Our theory suggests that un- 
conventional disorder effects may be responsible for the 
surprising absence of directed percolation scaling in at 
least some of the experiments. 
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